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Abstract: We extend previous computations of the second order QCD corrections to 
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state is massive. This allows accurate description of b — > ctv t decays. We review techniques 
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1. Introduction 

Inclusive semileptonic decays of S-mesons into charmed final states B — > X c Zi/; are bench- 
mark processes in £>-physics. These processes were studied extensively at ^-factories, LEP 
and the Tevatron ]l], |2], |J, [5], When lepton in the final state is electron or muon, most 
of the experimental data come from BABAR and BELLE, while measurements of inclusive 
rate for b — > X c tT> t transition are due to ALEPH and OPAL [@, g. At ^-factories only 
preliminary results on exclusive decay B — > Dtu t were reported Q. 

Theoretically, inclusive semileptonic decays of l?-mesons are well-understood, thanks 



to the Operator Product Expansion (OPE) in inverse powers of the 6-quark mass [10]. 
When results of experimental measurements for B — > X C W[ are compared with theoretical 
predictions, one is able to determine with high precision the bottom and charm quark 
masses, the CKM matrix element \V c b\ and the non-perturbative parameters of the heavy 
quark expansion [11, 12, 13]. Measurements of r lepton branching fractions by ALEPH and 



OPAL are used to constrain possible contributions of a charged Higgs boson to semileptonic 
-B-decays. 

An important issue in physics of semileptonic B-decays, intimately related to a very 
high precision achieved in experimental measurements, and the intention to utilize this 
precision fully, is the necessity to understand and control perturbative QCD corrections in 
heavy quark decays. There are two aspects of this problem. First, one needs to understand 



- 1 - 



the general structure of these corrections, to avoid large perturbative effects. A by- now 
standard example of this type is the recognition of a special role that short-distance low- 
scale quark masses play in avoiding large perturbative effects in quark decays [14|. Second, 
it is very important to provide explicit computations of QCD radiative corrections to 
quantities of direct relevance for experimental analysis. Large number of experimental 
results comes in the form of moments of lepton energy, hadron energy and hadron invariant 
mass with variable cut on the lepton energy. These are sufficiently complicated observables 
to make analytic computations, in particular in higher orders, impractical. On the contrary, 
if QCD effects are computed numerically for fully differential decay rates, any restriction 
on final state particles can be imposed. 

It is interesting to point out in this regard that, while one-loop corrections to the 
decay rate and a number of basic differential distributions for b — ► X C W[ decays were 
computed long ago |[^, fully differential decay rate through 0(a s ) was obtained only 
in 2004 [17, [l^, [l9|]. Early estimates of two-loop QCD corrections to the total rate for 
b — > X c eu e decays were given in Refs. |20|, |22|. Fully differential decay rate through 
0{a 2 s ) for b — > X c ev e was computed very recently in Ref. [23]. That calculation is based 
on the techniques developed in Refs. [24] for the computation of next-to-next-to-leading 
order (NNLO) QCD corrections to a number of processes in hadron collider physics. Those 
techniques were first applied to describe weak decays of charged fermions in Ref. [25], where 
two-loop QED corrections to the electron energy spectrum in muon decay were studied. We 
note that analytic results for semileptonic decay rate b — > X c lD[ and moments of the lepton 
energy without lepton energy cut were computed through 0(a 2 s ) in the form of m c /m{, 
expansion in Refs. p6| , p^ [. Analytic results reported in those references are very handy - 
especially, when compared to numerical calculations reported in Ref. [^3| and in the current 
paper - but it is not clear a'priori if they can be used for moments with lepton energy 
cuts, because the precision achieved in experimental measurements is rather high. Finally, 
we note that efforts are underway to compute 0(a s ) corrections to Wilson coefficients of 
leading non-perturbative operators, that contribute to moments in b — > X c Wi transitions 

There are two goals that we pursue in this paper. First, we extend the calculation 



reported in Ref. |23j by considering massive lepton in the final state. Once this is accom- 
plished, we have a fully differential description of b — > X c lv\ transition, where I can be 
electron, muon or r. Second, we present a much more detailed study of 0(a 2 ) corrections 
in semileptonic -B-decays, than what was published in Ref. j23|. 

We point out that since our computations are numerical, the presence of a massive 
lepton in the final state is a relatively minor complication, so that the extension of the 
calculation reported in [23| to b — » X c tv t is straightforward. This feature is particular to 
numeric computations since for analytic calculations any additional massive particle in the 
final state is a serious complication. 

We use our results for second order QCD corrections to b — > X c tv t transition to 
show that perturbative QCD corrections to the ratio B(b — * X c tv t ) / B{b — * X c lv{) are 
very small. This ratio was measured at LEP by the ALEPH and OPAL collaborations 
with decent precision. Non-perturbative corrections to this ratio are computed within the 
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framework of the operator product expansion in the inverse mass of the 6-quark |]29| , 30, 31] 
and are also rather modest. Hence, the ratio of the branching fractions provides a simple 
constraint on the values of bottom and charm quark masses. More generally, it appears 
that the ratio of two semileptonic branchings is an observable that can be predicted very 
accurately in QCD so that a good measurement of b — > X c tv t decay branching fraction at 
future l?-factories can be very interesting. 

We also provide a detailed study of the second order QCD corrections to b — > X c Wi 
for massless lepton. We compute the so-called non-BLM corrections to a large number 
of moments of different kinematic variables, for various charm and bottom quark masses 
and for different values of the lepton energy cut. These results can be used to create 
interpolating functions for 0(a 2 ) non-BLM corrections, for the use in fits to moments in 
semileptonic B decays. 

The remainder of the paper is organized as follows. In Section ||, we describe technical 
details of the computation. We begin by explaining the phase-space parametrization that 
we employ and then discuss details pertinent to our computation of two-loop virtual, real- 
virtual and double-real emission corrections. In Section |3| we present the results of the 
calculation. We tabulate large number of non-BLM corrections to various moments for b — > 
X c eD e decay and discuss their potential impact on the results of global fits of semileptonic 
-B-decays. Then, we present 0{a 2 s ) QCD corrections to the decay rate b — > X c tv t . We 
conclude in Section [|. 



2. Technical details 



In this Section we discuss technical aspects of the computation of the QCD radiative 
corrections to b — > X c Wi, where lepton can be electron, muon or tau. We would like to 
develop fully numerical approach to the computation of these radiative corrections, since 
this is the only known way to achieve flexibility, required for the description of 6-decays. 
Let us stress that such numerical computations would have been very straightforward if 
not for divergences that occur both in virtual loops and in real emission processes, when 
radiated gluons or massless quarks become soft or collinear to other particles. 

The presence of singularities makes it necessary to develop techniques to extract them, 
before proceeding to the numerical computations. This is accomplished with the help 
of the sector decomposition [p^| . Sector decomposition can be applied to integrals with 
complicated polynomials in denominators, to find changes of variables that factorize all 
singularities of the integrand. Given sufficiently complicated polynomials, such variable 
transformations can not be established globally and one needs to split the integral into 
many "sectors", where such variable transformations can be accomplished. It should be 
stressed that the sector decomposition technique is algorithmic and, hence, can be easily 
programmed - one does not need to examine all the integrals that appear in the problem 
to find suitable changes of variables. 

Calculation of any physical quantity requires that squares of matrix elements are in- 
tegrated over phase-space, allowed for final state particles. It is important, to choose the 
phase-space parametrization which is sufficiently simple to avoid proliferation of terms in 
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the process of sector decomposition. We therefore start with the detailed discussion of how 
the phase-space can be parametrized. After parametrization of the phase-space is fixed, 
we explain how various parts of the NNLO computation are performed. 

2.1 Phase-space parametrization 

We need to consider phase-space parametrization for the following processes: i) Born 
b — ► cWf, it) single-gluon emission b — > cWig; Hi) double-gluon emission b — > cWigg. In 
the latter case, the parametrization is also valid for quark emission processes b — * cWiqq, 
for massless quarks. We do not consider final states with three charm quarks since this 
contribution is suppressed, for realistic bottom and charm quark masses. 

2.1.1 Born process: b — > clv\ 

We begin with the phase-space parametrization for Born process b — * clvi. We label 
momenta of particles and their flavors in the same way so that, for example, b refers to p&, 
the momentum of the 6-quark, where appropriate. The differential phase-space for final 
state particles reads 

dLips LO = \dc\[dl\\dv\{2it) d ^ d) (b-c-l-Pt), (2.1) 

where the integration is performed in d-dimensional space, with d = 4 — 2e. The integration 
measure is defined as 

m = , d -L ■ ( 2 - 2 ) 
(2tt) 2p 

The phase-space decomposition suitable for a three-body decay is carried out by assuming 
a sequence of two two-body decays. First, the b quark decays into an off-shell VF-boson 
and the charm quark, then the virtual H^-boson decays into the lepton I and the neutrino. 
We write 

dLips LO = — [dc][dW](2TT) d 6^ {b-c-W) [dl][dui](2Tr) d 6^ (W-l-ui). (2.3) 
2ir 



We now parametrize all entries in Eq. (|2.3D in such a way that the integration region 
is a unit hypercube. We begin with the leptonic phase-space. We note that the leptonic 
phase-space is universal for all processes involved in the NNLO computation. Hence, we 
can choose to calculate it in four dimensions and neglect all its dependence on e. We 
therefore write 



VV fib 

dLips^^ = o _ TTr , 1 dA 7 dA 8 , 0<A 7)8 <1. (2.4) 



W' 2 - m 2 
8ttW 2 

The physical meaning of the two parameters Aj^ is as follows - they describe the polar angle 
cos 9i = — 1 + 2A7 and the azimuthal angle 4>i = 2tt\% , that parametrize lepton momentum 
in the rest frame of the VF-boson, relative to the direction of the W momentum in the rest 
frame of the decaying 6-quark. 

The remaining phase-space that describes b — > c + W transition can be written as 



dW 2 „ Um b -m c ) -m 2 )p d c 3 ^ 

— dlapw^ = — dx 1; (2.5) 
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where dO c is the (d— l)-dimensional solid angle that describes direction of the charm quark 
momentum p c in the 6-quark rest frame. The remaining kinematic variables in the 6-quark 
rest frame are written as 



max / , , ">r '"l 



E c = m c + {Ef 3 * - m c ) (1 - x x ), S c max = 1 + 



2 V ' ml 



E w =m b -E c , P w = ^E^-W 2 , W 2 = m 2 + m 2 c - 2m b E c . (2.6) 

It is straightforward to write the energy of the lepton and the angle between the lepton 
and the W in the 6-quark rest frame 

El _ % ( 1 + £ V *Z U - ^ (1 - 2A 7 ) , cos 9iw = 2E " E <- W2 ~ m - . (2.7) 



2 V W 2 J 2 V W 2 J y " 2P wP i 



In Eq.( |2.7| ), pi = y Ef — mf denotes the lepton three- momentum. 

The above formulas provide sufficient information to compute scalar products of the 
four-momenta of all particles that appear in the leading order calculation. Indeed, scalar 
products that involve the decaying b quark are straightforwardly expressed in terms of 
particle energies in the 6-quark rest frame. The relative angle between the direction of the 
charged lepton and the direction of the charm quark is easily related to Oiw- Neutrino 
four-momentum is given by v\ = W — I which implies that neutrino momentum does not 
lead to independent scalar products. The phase-space parametrization described in this 
subsection, is employed in all parts of the calculation where the three-particle phase-space 
enters. In NLO and NNLO computations, this occurs when contributions to the decay rate 
due to one- and two- loop virtual corrections, respectively, are calculated. 



2.1.2 Single gluon emission process: b — ► clvig 

Next, we discuss the phase-space parametrization for the process b — ► cWig, with four 
particles in the final state. When the energy of the emitted gluon becomes small, the 
corresponding matrix element diverges. Good parametrization of the four-particle phase- 
space should factor out the dependence on the gluon energy, so that extraction of infrared 
divergences occurs easily. We write 

dLips NLO = [dc] [dg] [dl] [din] (27r) d 5^ (b-c-l-^-g), (2.8) 

and decompose it into quark and lepton phase-spaces, by introducing the four-momentum 
of the W boson. This leads to 

dW 2 

dLips NLO = — r dLips 6 _ KS+tf+w -dLips W r_, +Si . (2.9) 

Parametrization of the lepton phase-space is the same as for the Born process, described 
earlier. The quark phase-space is different. To arrive at a suitable parametrization, it is 
convenient to integrate over the three-momentum of the H^-boson and then over the energy 
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of the gluon, to remove all delta-functions. We introduce three variables < £1,2,3 < 1 and 
write 

(777-6 — w-c) 2 — mf (1 — X\) X 2 

E c = m c + (Er x -m c ) Xl , E g = ^— — — ^ — , (2.10) 

2(m b - E c + p c cos 9 cg ) 

W 2 = TJlf + (7775 — 777 c ) 2 — mf (1 — X\) (1 — X 2 ) , COS0 C g = — 1 + 2^3- 

All other angles can be derived. For example, the angle between the gluon and the W in 
the 6-quark rest frame reads 

-Eg -p c cos9 cg 
cos6 gW = y —— y -. (2.11) 

fy 



The only other angle that we need in order to fix all independent scalar products is the 
angle between the lepton and the gluon, in the 6-quark rest frame. We obtain 

cosOig = cos 8gw cos 9 iw + sm8 g w sin Oiw cos <pi, (2-12) 



where cosfyw is given in Eq.(2.7) and (f>i = 2-7rA8. Finally, we obtain the following 



parametrization of the NLO phase space 

3-2e 



dLips NLO Q 2 T(2-2e) Pc 2e {{m b - m c ) 2 - mf 



ft dx,dA 7 dA 8 28 ~ 267r ^f {d ~ l) r2 (! " £ ) m * ^ - E c + Pc cos 9 cg f~' 1 - 
i=l 



x (l - (1 - *!) 2 - 2e (1 " • (2-13) 

We use dLips^Lo m NLO calculations, as well as for dealing with real-virtual corrections 
in NNLO calculations. 

2.1.3 Double gluon emission process: b — > clv\gg 

Finally, we discuss the parametrization of the five particle phase-space that is needed for 
the description of the double real-emission processes, such as b — ► clvxg\gi or b — > cWiqq. 
We introduce six variables Xj = i..g, that satisfy < X{ < 1 and use £1,2,3 to parametrize 
energies of the charm quark and of the gluons, and the W invariant mass 

E c = m c + (£ c max " m ) (1 - xi), W 2 = mf + {(m b - m c f - m\)x x {\ - x 2 ), (2.14) 

((777 6 - 777 c ) 2 - mf)xiX 2 X 3 ( (m b - 777 c ) 2 - 777 2 ) X X X 2 (1 ~ X 3 ) 

E n , — — ; — — ; , E, 



" 1 2(mfe — E c + p c cos 9i c ) ' 92 2(7775 — E c + p c cos 9 2c ) — E gi (1 — COS I 



12 



We use £4,5,6 to parametrize the relevant angles. We note that, in order to handle collinear 
singularities related to g* — > gg splitting, it is useful to have a simple parametrization of 
the relative angle between the two gluons. Therefore, we choose the z-axis to be aligned 
with the momentum of the gluon g\ ; we choose the x-axis in such a way that the gluon g 2 
is in the z — x plane. This fixes the global reference frame. Then, we introduce the relative 
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angles between two gluons, gluon g\ and the charm quark and the azimuthal angle of the 
charm quark 4>c 

cos #12 = 1 — 2x4; cosmic = 1 — 2x5; cos 4>c = 1 — 2 sin 2 (J~^J ■ (2-15) 

Given those angles, we can find all other angles between momenta of different particles. 
For example, the angle between the charm quark and the gluon 52 reads 

cos 9 C 2 = sin 612 sin 9\ c cos <j> c + cos 9\2 cos 0\ c . (2-16) 

The angle between the M^-boson and any other particle is computed from momentum 
conservation W = b — c — g% — Q2- For example, the angle between the VF-boson and the 
gluon g\ is given by 

co S 9 lw =- E ^- E ^Z dl2 ~ PcCOS9lC - (2-17) 

Pw 

The relative angle between the charged lepton and any other particle is derived in a similar 
way. For example, the angle between the gluon g\ and the lepton I reads 

cos On = cos Oiw cos Q\W + sin#ivp sin#^ cos(^). (2-18) 

Finally, we express the phase-space for b — ► clvig\g2 decay through the appropriate 
variables and obtain 



dLips NNLO ft 3 r 2 (2 - 2e)r(l - 2e) 



,2 



1 " T77^ (2-19) 



6 2 n - 4e (27r) M - 3 r 4 (l-e)r 2 (l/2-e) V W 2 

]1 dxidA 7 dA 8 

i=l 

((mft — m c ) 2 — mffi~^' L, p)r 2e 
m b (m b - E c + p c cos 6 lc ) 2 - 2t (m b - E c +p c cos 6 2c - E gi (l - cos 6>i 2 )) 2 ~ 2e 
xx\-^xl-^x\- 2 \l - x 3 ) 1_2e xr e (l - x 4 )- £ x 5 - e (l - x 5 )- e x 6 - e (l - x 6 )- e . 

Having discussed parametrization of phase-spaces that we employ in the NNLO com- 
putation, we continue with the description of technical details relevant for the computation 
of 0{a 2 ) QCD corrections to b — ► clui transition. There are three distinct components that 
need to be addressed - two- loop virtual corrections, virtual corrections to single real emis- 
sion process and double real emission corrections. Since these components require different 
techniques, we describe the relevant details in the following subsections. 

2.2 Two-loop virtual corrections 

We begin with the discussion of how two-loop virtual corrections are computed. There are 
twelve two-loop diagrams; examples are shown in Fig.[l]. These diagrams are complicated 
because they involve several mass scales, m 2 ,m 2 ,W 2 ,m 2 as well as complicated tensor 
integrals, e.g. due to spin correlations of final state leptons with bottom and charm quarks. 
These features make analytic computations impractical. However, we can compute those 



diagrams numerically using the method of sector decomposition [32|. We point out that 
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c) d) e) 

Figure 1: Sample two-loop diagrams that contribute to b — > c + W transition. 



application of sector decomposition is simplified in 6-decays, since the two-loop diagrams 
do not develop genuine imaginary parts. In principle, sector decomposition method was 



extended recently |33| 34 1 to deal with problems where imaginary parts do appear but it 
is a welcome feature of the problem at hand, that we do not need to deal with additional 
complications. 

Hence, the primary issue in the calculation of two-loop Feynman diagrams that we have 
to address is the efficient choice of Feynman parameters, to reduce the amount of sectors 
that are created in the process of sector decomposition. It is also important to perform 
integration over loop momenta in such a way that computation of tensor integrals does not 
introduce kinematic singularities. It turns out that a very simple and fairly efficient way to 
deal with tensor integrals is to integrate over two loop momenta sequentially. To illustrate 
this procedure, we consider the planar two-loop diagram, shown in Fig. la. This diagram 
can be represented as a linear combination of tensor integrals 



d d h r d d k 2 {fciMMj 

(27r) d J (2vr) d D^D^DM 



where k\ and k 2 are the loop momenta, {q}i = q^.-.q^ is the rank-i tensor, composed of 
the relevant loop momenta and -D m =i..6 are inverse Feynman propagators that appear in 
the planar diagram. They read 

D 1= kf, D 2 = kf + 2kib, D 3 = kf + 2k lC , (2.21) 
D 4 = kl D 5 = k 2 12 + 2k 12 b, D 6 = k 2 12 + 2k 12 c, (2.22) 

where k% 2 = k\ + k 2 is used. To integrate over momentum k 2 , we introduce two Feynman 
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parameters and write 



A - -s/ldW ^y-.?', . (2.23) 



(2vr) d D 4 D 5 D 6 P Ai J (2vr) d ( K 2 - A(h)f ' 

where [dx 2 dx 3 ] = dx 2 dx 3 #(l — x 2 — x 3 ) is the integration measure and 

K = k 2 + Q, Q = (k 1 + b)x 2 + (k 1 + c)x 3 , (2.24) 
A(fci) = Q 2 - (/c 2 + 2ife a 6) x 2 - (/c 2 + 2fc]c) x 3 . (2.25) 

Integration over the shifted loop momentum K is standard and can be easily performed 
for arbitrary rank tensor. The important point is that the higher the rank of the tensor, 
the smaller the power of the function 1/A is in the resultant integral. Since numerical 
integration is mostly problematic because of infrared divergences, we should be looking at 
the most infrared-singular integral, which is provided by the JC-less term in the numerator 
of Eq.( 2.23| ). For such integrals we find 



2 



1 

/\, 1 f d d K W}j »r(l + e) /\, , . \Q}j 
/ dx 2 dx 3 / 7 — -7 o = , ZJ / dx 2 dx 3 — 1 , ■ ■ (2.26 

Note that integrals with higher powers of K in Eq.( |2.23| ), can, after K- integration, be 



written as a linear combination of the integrands in the right hand side of Eq.(2.26), by 
multiplying and dividing the integrand by appropriate powers of A(fci). The important 
point is that additional powers of K in Eq.( 2.23|) do not generate yet higher powers of 1/A 



in Eq.(B) 



The next step requires integration over fci; to do that, it is convenient to change 
variables x 2 = A1A2, x 3 = Ai(l — A 2 ). The integral over k\ becomes 

f d\ 1 d\ 2 f {k} tl 
13 11 ~ 1 Af(l - Ax)i+« J (2nY Dl D 2 D 3 A^ ' 1 ] 



where 



and 



A = A; 2 + 2^2--^-, (2.28) 
1 — Ai 



Q 2 = A 2 6+(1-A 2 )c. (2.29) 



Since A is a polynomial in k\, the integration over k\ can be performed in, essen- 
tially, the same way as what was described above in regard with k 2 integration. Integrals 
with strongest infrared divergences are the ones without additional powers of k\ in the 
numerator. The corresponding scalar integral reads 

r(2 + 2 £ ) f" Ap(l - Ai) 1+£ A 4 A|(1 - A5) 2 
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where 



F = (1 - Ai)(Q 3 A 4 (l - A 5 ) + Q 2 A 5 ) 2 + Q^AiAg 



(2.31) 



and 



Q 3 = A 3 6+(1-A 3 )c. 



(2.32) 



It is clear from Eqs.( |2,30 , 2.31| ) that the function \/F 2+2<L develops overlapping sin- 
gularities at the integration boundaries; for example F = for Ai = l,As = and for 
A4 = 0, A5 = 0. To disentangle those singularities, we employ the technique of sector 



decomposition 32]. To this end, we map all the singularities to the origin by splitting 
the integration region into two intervals [0,1/2] and [1/2,1], for each Aj, and then change 
variables Aj —* Aj/2 and Aj — > 1 — Aj/2 in the first and second interval, respectively. The 
sector decomposition is then applied to the integrand; this allows us to find a sequence 
of variable transformations that factorize all singularities. Once singularities are factored 
out, for each tensor integral we get expressions of the following type 



a € sect 



Iij = ^2 dAidA 2 ...dA n 



A^(A!,A 2 ...A n ) 



UXl+ a ^D9(\ 1 ,X 2 ,...X n ) 

1=1 



(2.33) 



where all functions -A^"({Aj}) and Dfj({Xi}) are finite throughout the integration region. 
Hence, all the singularities of the integrand are in explicitly factorized form and it is easy 
to obtain integrable expressions by employing the plus-distribution prescription 



where 



f 

Jo 



\l+c 



dXf (A) 



—6(X) + Y] 
ae ^ 



n=0 



ln n (A) 



(-ae) 



"ln n (A)~ 


=/ 


A 


+ Jo 



/(o) 



In™ (A) 



(2.34) 



(2.35) 



We are now in a position to sketch all the steps that we go through to carry out a calcu- 
lation of a planar two-loop diagram. After the planar two-loop diagram is multiplied with 
the complex conjugate tree- level amplitude and summation over polarizations of external 
particle is performed, we use Form |3j| to integrate over the loop momenta, following the 
procedure that we just described. As explained earlier, we can always write the result in a 
form similar to that of the scalar two-loop integral provided that we allow for a polynomial 
function of Feynman parameters in the numerator. Since those numerator functions are 
finite, we do not need to know their explicit form and can treat them as generic finite func- 
tions in the process of sector decomposition. The sector decomposition procedure is coded 
up in Maple |^6|. After sector decomposition is completed, Fortran files that contain finite 
functions to be integrated and all the changes of variables that the sector decomposition 
procedure found necessary to apply to the integrand to factor out potential singularities, 
are automatically written out. 
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Computation of the non-planar diagram and the diagram with the three-gluon vertex 
is very similar to what is described above; all that changes is the Feynman parametriza- 
tion. However, the procedure has to be modified for diagrams with vacuum polarization 
insertions on a gluon line and for diagrams with self-energy insertions on bottom and 
charm quark lines. We begin with the discussion of the vacuum polarization diagrams with 
massless particles, e.g. gluons, quarks and ghosts. Such diagrams read 



' (2n) d J (2n) d kl{k l + k 2 f(kfj 1 {kl + 2k l h){k\ + 2k l c) 



An obvious issue here is the presence of two identical gluon propagators k^ 2 . As the 
result, in the limit k\ — > 0, the denominator in Eq, (|2.36| ) develops cubic, rather than 
linear, singularity. In principle, even in this case, it is possible to proceed along the lines 
described above for the planar diagram all the way through the application of the sector 
decomposition and factorization of singularities. However, the complication occurs in the 



process of the extraction of singularities using plus-distributions Eq.( 2.34 ), since a term 
that scales like x~ l ~ n for x — > 0, leads to an expansion that involves n-th derivative of a 8- 
function rather than the ^-function itself. As it turns out, this complication is unnecessary 
since one can analytically integrate over k 2 in any massless vacuum polarization diagram 
and observe that 

d d k 2 N(kx,k 2 ) 



(2w) d k\ (h + k 2 y 



kfNih), (2.37) 



thanks to gauge-invar iance. When dealing with massless vacuum polarization diagrams, we 
indeed integrate over k 2 analytically, cancel one of the \jk\ propagators and then perform 
numerical integration over k\. 

Clearly, a similar problem occurs also in the case of vacuum polarization corrections 
with massive quarks. In that case, however, it is harder to explicitly factor out the depen- 
dence on the loop momentum k±, to cancel the cubic divergence at small k\. For vacuum 
polarizations with massive quarks we adopt a different strategy - we subtract those vac- 
uum polarization loops at zero momentum transfer and use dispersion representation to 
connect the two-loop diagram with the massive fermion loop to a one-loop diagram with 
the massive gluon [37]. 

Non-integrable singularities appear also in diagrams with the self-energy insertion on 
the massive (b or c) lines; in this case they are caused by the square of the massive propa- 
gator which becomes nearly on-shell for small momentum of the virtual gluon. In variance 
with the case of the massless vacuum polarization, it is not possible to perform analytic 
integration over the loop momentum of the "self-energy" loop. To get around this problem, 
we use a particular integral representation for the quark self-energy diagram. We consider 
a self-energy diagram 



d d k 2 ^(p + h + mhii 
(27r) d ((p + k 2 y-mf)kr 



(2.38) 
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where mi stands for mj, or m c . We combine two denominators using Feynman parameters, 
integrate over the loop momentum and obtain 



i 



»r(e) f 2 2x-e l^iP + mi -px)^ 



(4vr)^ J dXX ~ e(P - m ^ \l- p 2 x/{p 2_ m 2 )r ( 2 - 39 ) 



This integral can be written through hyper geometric functions. To this end, we introduce 
two Dirac structures 

Ni =^(p + m i )j tM = (2-d)p + dm, N 2 = 7^ = (2 - d)p, (2.40) 

and write 



(p 2 - mf)' 



^ ^21 U, 1 - e; 2 - e, 



( 4 ^)d/2 



r(i-e) r / , P 



2 



r(2-e)^ zl V"^ ^p 2 -m 2 



r(3 — e) \ ' ' P 2 — mf 



(2.41) 



These hyper geometric functions are not suitable for the purpose of subsequent integration 
because the on-shell limit p 2 — > mf is at infinity. To take care of that, it is useful to employ 
an identity that relates the hyper geometric functions with the argument z and z/(z — 1). 
If we perform this transformation and go back to the integral representation, both of the 
hypergeometric functions (HGFs) become 

1 

dx 9 [ 7, , (2.42) 

jr — mf I x 





- - m 2 



where a, b are some integers. It is now straightforward to subtract from and add to the 
integrand in the right hand side of Eq.( 2.42j ) its value at p 

s -ar 1 _ x \b-2e p 2_ m f r x £ - a+1 (l - x) b ~ 

dx =i / = g-t- / dx ^ / mr 2 B{e -a + 2,b- 2e). 

p A — mf/x m\ J p^ — mf/x 







(2.43) 

The first term in the right hand side of Eq.( 2.43| ) can be inserted in a two- loop diagram 
since the pre- factor p 2 — mf cancels one of the problematic massive fermion propagators. 
The second term is a constant; it can be combined with the mass counter-term contribution 
to cancel quadratic singularity and, after that, calculated explicitly. 

2.3 Mixed real-virtual corrections 

In this subsection we discuss the calculation of one-loop radiative corrections to single gluon 
real emission amplitudes. At first sight, these corrections may look much simpler than the 
two-loop virtual corrections discussed previously, since they involve only one-loop virtual 
diagrams and the final state with relatively low multiplicity. It would appear therefore that, 
technically, they fall into a category of well-established next-to-leading (NLO) calculations 




a) b) 



Figure 2: Sample diagrams that describe one-loop corrections to b^c + W + g transition. 



pql - Unfortunately, this is not quite true since, in contrast to standard NLO computations, 
the emitted gluon in the final state can become soft, invalidating applicability of NLO 
computational techniques. Therefore, real-virtual corrections require careful study. 

There are two strategies that one can pursue to deal with the real-virtual corrections. 
One can use Passarino-Veltman tensor reduction technique and integration-by-parts 
identities j^]] to reduce real-virtual corrections to one-loop scalar integrals. Then, one 
can attempt to extract singularities that appear when the energy of the gluon in the final 
state becomes small. While this approach was used in a number of calculations [^4|, |25|] , it 
rapidly becomes impractical with the increase in the number of particles in the final state. 

A flexible method should be based on numerical computations and it seems that Feyn- 
man parametrization of one-loop virtual corrections and subsequent application of sector 
decomposition to both Feynman parameters and the energy of the emitted gluon is a 
straightforward thing to do. The only problem with that approach is that one- loop cor- 
rections to real gluon emission do develop imaginary parts, even when all parameters in 
the integrand are real. Technically this happens because of the singularity on the real 
integration axis which is regulated by the +i0 prescription. While it is easy to implement 
such a prescription in analytical computations, it is difficult to do so in a fully numerical 
approach. 

It turns out that there is a simple way to avoid the issue of the imaginary part in this 
problem. To this end, we observe that, for any Feynman diagram of the real-virtual type 
that contributes to b — > c transition, it is possible to choose integration variables in such a 
way that the integration over at least one variable is of the form 

1 

dx / dyF(x,y), F(x,y) = ^rr- (2.44) 

In Eq. fl2.44j ) a: is a collection of other variables involved in the computation of the integral, 
A and B are some functions of those variables which satisfy B({x}) > A({a;}) for all x 
and nj = i..4 are integers. We do not have a proof that such a parametrization is possible 
for real-virtual corrections under all possible circumstances, but we find empirically that 
it exists for b — > c transitions. 
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The problem with the numerical evaluation of the integral in Eq.( 2.44 ) is that, for 
rt3 > 0, it becomes singular at y = A(x) / B{x) < 1, so that this singularity occurs in 
the middle of the integration region. Such singularity can not be included into the sector 
decomposition framework in a straightforward way. To deal with this problem, we rewrite 
Eq.( |2.44} ) in the following manner 



dx / dyF(x,y)- dx dy F(x,y) 



(2.45) 



and observe that the first integral can be computed analytically, while denominator of the 
function F(x,y) in the second term is sign-definite. Changing variables y — * 1/y in the 
second term in Eq.(p,45]), we obtain the integral that is amenable to sector decomposition. 



We now illustrate this general discussion by considering explicit examples. To set the 
stage, we begin with a simple case, where the imaginary part problem does not occur. 
This happens for all diagrams where the gluon is emitted from the 6-quark line. For our 
example, we consider diagram Fig. ^a. Interference of this diagram with the tree amplitude 
that describes radiative decay of the 6-quark b — > c+W+g, contributes to 0{o? s ) correction 
to the decay rate. We consider a scalar integral associated with this loop diagram 



1 



d d k 



1 



T(l + e) J in d / 2 k 2 (k 2 + 2k{b-g)- 2bg) (k 2 + 2kc) ' 
We introduce Feynman parameters and integrate over the loop momentum to obtain 



(2.46) 



I 



[dx 2 dx 3 



((6 - g) x 2 + cx 3 ) 2 + 2bgx 2 



1+6 



/ 



dAidA2 



a; 



[A(Ax,A 2 )] 



l+e' 



where 



A(Ai, A 2 ) = Ai ((b - g) A 2 + c(l - A 2 )) 2 + 2bg\ 2 . 
Writing b — g = c + W, we find 

A = Ai(m 2 + 2cW\ 2 + W 2 \ 2 2 ) + 2bg\ 2 > 0, 



(2.47) 

(2.48) 
(2.49) 



for all values < Ai, A 2 < 1. This representation for A(Ai, A 2 ) is instructive since it shows 
how new overlapping singularities appear when the emitted gluon becomes soft. Indeed, 
for non- vanishing gluon energy, A(Ai, A 2 ) vanishes for Ai = A 2 = 0. On the other hand, if 
the gluon is soft g — » 0, A(Ai, A 2 ) vanishes for Ai = and any value of A 2 . To take care of 
all possible cases, we employ explicit parametrization of the gluon energy in the 6-quark 
rest frame, as described in the previous Section, and perform sector decomposition of the 
amplitude squared treating Ai,A 2 and bg on equal footing. This allows us to extract all 
singularities associated with vanishing of both, the loop momentum and the momentum of 
the gluon in the final state. 

We now turn to the description of a more difficult case which occurs when, in the 
one- loop amplitude, the gluon is emitted from the charm quark line. A representative 
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diagram is shown in Fig.^. All the problems that appear in this case can be illustrated by 
considering the scalar integral as an example. We have 



1 f d d k 



r(l + e) J iir d / 2 k 2 {k 2 + 2kb) {{k + c + g) 2 - m 2 ) 
[dx 2 dx 3 ] 



2 1 1+e ' 

(bx 2 + (g + c)x 3 ) - 2gcx 3 



(2.50) 



Changing variables x 2 = Ai(l — A2) and X3 = A1A2, we obtain 

dAidA 2 ^ r— = I dAidA 2 ^ t— , (2.51) 

2 (AA!-2c 5 A 2 ) 1+e J 2 A^(Xt-0 1+e 



where 



A = (b-W\ 2 )\ e = ^<l. (2.52) 



It is clear from Eq.( |2.5l| ) that there is a singularity at Ai = £, i.e. in the middle of 
the integration region, which can not be dealt with using the sector decomposition. To 
transform Eq.(2.51) into a form suitable for numerical integration, we proceed along the 



lines described in the beginning of this Section. Consider integration over Ai. Allowing for 
a general form of the integrand, we re-write 



1 00 

a— e 



d V-V< = J d V-6^ " J dAl (A!-0^ 

1 

= e+a - b - 2er{1 _ b _ e) L 1)6+£ r(a + l-e) + T(2e - a + b - 1} 



T(2 + a-6-2e) r(e-a) 







"^ (i-^)^ - (2 ' 53) 



The remaining integrand is sign-definite, because £ < 1. The problem of the singularity 
in the middle of the integration region has been handled by an analytic integration; the 
remnant of this problem is the imaginary part in (— l) b+e . It is now very straightforward 
to apply the sector decomposition to the remaining integral in Eq.( 2.53| ), to extract the 
infrared singularities. 

The above technique is applicable to all real-virtual diagrams. The first step - finding 
Feynman parametrization that is linear in (at least) one variable - requires careful inspec- 
tion of each diagram individually but, once such parametrization is found, remaining steps 
are easily accomplished using algebraic manipulation programs. A slightly different ap- 
proach is required for real-virtual contributions related to self-energy insertions on massive 
fermion lines where quadratic singularities appear. We deal with those singularities using 
parametrization of self-energy diagrams described earlier in this Section in the context of 
two- loop virtual corrections. 
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w* w* 

Figure 3: Diagrams that describe double real emission process with a massless qq pair in the final 
state. 



2.4 Double real emission corrections 

Finally, we address the computation of the double real emission corrections. Because we 
deal primarily with massive quarks, the majority of diagrams develop only infrared singu- 
larities. We can easily extract those using explicit parametrization of gluon energies in the 
6-quark rest frame, discussed earlier in this Section. The only exception to the "no collinear 
singularity" rule comes from diagrams where an off-shell gluon splits into two gluons or 
a qq pair. In principle, since we use the relative angle between two gluons (or massless 
quarks) as one of the primary variables for the phase-space parametrization, extraction of 
(potential) collinear singularity is also straightforward. However, the complication arises 
because of the necessity to deal with the extraction of quadratic singularities in some of 
the diagrams. The problem and the solution is best illustrated by considering diagrams 
with a massless qq pair in the final state. 

Consider two diagrams that contribute to the process b — ► clvi + qq, shown in Fig. 3. 
Upon squaring these diagrams, we find that the gluon splitting g* — > qq leads to a structure 
that involves square of the gluon propagator 

(qq) 2 • ( } 

An obvious problem with this result is that the singularity associated with the collinear 
limit is power-like; hence, we can not disregard the numerator in the process of the sector- 
decomposition since it provides the necessary qq scalar product to soften the collinear 
singularity. The question is how to parametrize the momenta q, q to enable easy extraction 



of the scalar product qq from the numerator in Eq.(2.54). To deal with this problem, 
we need to exploit the fact that, when collinear singularity occurs in those diagrams, the 
momenta of q and q become parallel to each other. We use this feature to parametrize the 
scalar products of q and q momenta with, say, the charm quark momentum in the following 

way 

cq = ^-cq + Jxlxxx^cs- (2.55) 
bq 

Here X4 describes the angle between q and q, cos #12 = 1 — 2x4. The important point is the 
factor yjxl in the second term on the right hand side in Eq.( [2.55 ); it is crucial for regulating 
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collinear singularity. On the other hand, the (complicated) function A C g does not need to 
be made explicit in the matrix element to cancel collinear singularity. Nevertheless, we 
give it here for completeness 

^ _ _ 4p c {(m b - m c ) 2 - mf) (1 - x 3 ) (^/x4Cos#i c - y/1 - x 4 sin 9\ c cos <j) c ) ^ 
2(m b - E c + p c cos 02c) - E 9l (1 - cos 9 12) 

Very similar manipulations are needed for scalar products that involve the momentum of 
the charged lepton and the momenta of q and q. 

3. Results 

In this Section, we discuss the two-loop QCD radiative corrections to b — > X c li>i transitions 
and present results for a large number of moments, relevant for the experimental analysis 
and semileptonic fits. We begin by writing the decay rate b — > X c li>i as 

where d-Fj stands for the differential decay rate, at leading, next-to-leading and next-to- 
next-to-leading order, respectively. The strong coupling constant a s = a s (m&) is defined 
in the MS scheme in the theory with three massless flavors; it is renormalized at the value 
of the 6-quark mass. The lepton and hadron moments are defined as 

L n (E cnt ) — , (3.2) 

„ ,„ , <(K ~ OKY (E h /m b y 9 (Ei - E cnt ) dT) 

\dl 0/ 

where Eih are the lepton and hadron energies in the 6-quark rest frame and is hadronic 
invariant mass. Also, 

^^Mj4 iFo , (3.4) 
192-k- 3 

and (..) implies integration of the corresponding quantity over available phase-space 1 . The 
calculation is performed in the pole mass scheme. For numerical integration, we use Vegas 



f40f| , implemented in the Cuba library [41]. We treat the axial current as suggested in 
Ref.gg. 

The lepton and hadron moments can be computed in an expansion in the strong 
coupling constant 



L„ LT + f L« + ) 2 ( PoLt BLM) + + (3.5) 



1 Available phase space at the parton level is determined by the value of the 6-quark mass. Cuts on the 
phase-space are shown explicitly in Eq. (J3. l[) . 
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Table 1: Non-BLM corrections to lepton moments L\ in dependence of r and £. Vegas integration 
errors are shown in brackets. 



where /3q = 11 — 2/3nj, and nj = 3 is the number of quark flavors that are treated as mass- 
less in the computation. Next-to- leading order and BLM corrections [43| to any kinematic 
distribution in b — » X c lT>i transition are known [15, 16, 17, 18, 19|. Non-BLM corrections 



(2) (2) 

Ln , H-j for massless lepton were computed very recently and their detailed investiga- 
tion is not available. In the remainder of this Section, we study non-BLM corrections to 
b — > X c eu e decay in detail. Then we describe results for the inclusive rate b — > X c tv t . 



3.1 Non-BLM corrections and moments of b — > X c ev e decays 

In this subsection, we study corrections to semileptonic decay b — > X c Wi, where I is the 
massless lepton, in dependence of the bottom and charm quark masses, and the lepton 
energy cut. Because lepton and hadron moments defined in Eq,(|3.5|) and Eq.( |3.6|) are 
dimensionless, they depend on the two ratios of the three dimensionfull parameters. We 
choose r = ruc/m^ and £ = 2E' cut /m£, as independent variables. In Table 1 we show non- 
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BLM corrections to lepton moments Lo,i,2,3 for a number of r and £ values. In Tables 2,3 
results for hadron moments are given. These results can, in principle, be used in global fits 
of semileptonic decays where b and c masses are parameters that need to be fitted. 

It was observed in [23] that second order QCD corrections to 6-decays do not depend 
strongly on kinematics and it is interesting to further explore this observation. To this 
end, we may conjecture that non-BLM corrections to moments are given by constant, 
^-independent renormalization factors of the leading order moments. If true, this renor- 
malization factor can be determined from Refs. ^7|, where lepton and hadron energy 
moments are analytically computed for zero lepton energy cut. We will construct the inter- 
polating function for lepton energy moments, following this observation. At leading order, 
the lepton energy moments are given by 



F o (r,0)' 



(3.7) 



where 



dx 



2x 2 (x r 



(i- x y 



(6 — 6x + xx m + 2x 2 — 3x m ) 



(3.8) 



In Eq. (|3.8[) x m = 1 — r and the integration is over x = 2Ei/mb- Since the integration here 
is elementary but the resulting formulas are lengthy, we do not present the results of the 
integration here. We introduce the interpolating function by defining 



L 



L\ ; (r,0) 



(3.9) 



so that the normalization of the non-BLM correction to the moment is fixed by its value 
at zero lepton energy cut and the shape is taken to coincide with the leading order shape. 
The interpolated moments L^' m (r, £) are given in Table 4. Comparing computed and 
interpolated moments, we observe that L^' m (r, £) provides excellent approximation to 

(2) 

L\ (r, £) for small values of £. However, the agreement becomes progressively worse for 
larger values of £. For example, a typical deviation between the interpolated and the 
explicitly computed non-BLM moments for £ = 0.7 and m c /mb = 0.28 can be as much as 
twenty percent. Finally, we point out that a very similar behavior is observed for non-BLM 
hadron energy moments. 

Note that by increasing the cut on the lepton energy, the phase-space is restricted 
to the region where soft gluon radiation becomes relatively more important and, hence, 
the dynamics of the final state changes with the increase of the cut on the lepton energy. 
It is therefore clear that for moments defined in Eqs.( |3.5| , |3.6[ ) perturbative corrections at 
different lepton energy cuts are not correlated, i.e. different physics becomes important 
for different values of the lepton energy cut. On the other hand, we point out that the 
moments measured in experimental analysis correspond to ratios of L-moments defined in 
as we explain below, this difference is essential for understanding importance 
of QCD radiative corrections in the global fits. 
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Table 2: Non-BLM corrections to hadron energy moments in dependence of r and £. Vegas 
integration errors are shown in brackets. 



We turn to the discussion of the potential impact that computed corrections may have 
on the extraction of fundamental quantities in heavy quark physics, such as \V c b\, nib, m c , \Xq 
etc. from global fits to semileptonic moments. We stress that we do not attempt to perform 
a fit to data on semileptonic moments, leaving this task to experts. However, we find it 
instructive to illustrate shifts that may be expected in the values of, e.g. the \V c b\ and the 
6-quark mass, if non-BLM corrections are included. 

We begin with the discussion of the CKM matrix element \V c b\- Since the |V^| is 
obtained from the normalization of the partial decay rate, it is mostly sensitive to QCD 
corrections to the moment Lq. The non-BLM corrections to that moment for £ = 
and various values of r are shown in Fig^. We see that, for realistic ratios of quark 
masses, the non-BLM corrections to Lq are between 2 and 1.5 percent. Since experimental 
measurement fixes iV^pLn, one can expect that \V c b\ changes by about —1 percent, when 
non-BLM corrections are included. This is compatible with the uncertainty in \V c b\ as 
currently estimated. 

BABAR collaboration measured a number of lepton energy moments for b — > clui 
transitions with high precision Q. For the illustration, we use their measurement of the 
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Table 3: Non-BLM corrections to hadron invariant mass moments in dependence of r and 
£. Vegas integration errors are shown in brackets. Integration errors are not shown if they are 
significantly smaller than one percent. 



moment Aii for two values of the lepton energy cut 

f Ei ar , 

u {T? Scut _ J 1437.6(4.0) (5.7) MeV, E cut = 0.6 GeV; 

Ml{E ^ ] ~ -J^T ~ { 1773.7(1.9)(1.1) MeV, E cut = 1.5 GeV. (3 - 10) 



C \ 1 1 



These (and other) results are used to extract the following values of the bottom and charm 
quark masses in the kinetic scheme [44] m?* = 4.55(5) GeV and mf = 1.08(7) GeV, where 



we combined all uncertainties in quadratures. In terms of lepton moments computed in 
this paper, and neglecting non-perturbative contributions, we find 

M x {E cnt ) = m b Mf, Mf = (3.11) 

The non-BLM corrections to lepton moments were not accounted for in || . To estimate 
their impact, we compute Lq and L\ for nib = 4.55(5) GeV and m c = 1.08(7) GeV and 
-Scut = 0.6 and 1.5 GeV. The results are given in Table 5. Expanding Eq.( |3.1l| ) around 
mi, = + Snij,, to account for the shift in the 6-quark mass, induced by including non- 
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Figure 4: Non-BLM correction to the total rate, in percent, for zero lepton energy cut, as a 
function of the charm quark mass to bottom quark mass ratio. We use a s = 0.22. 



BLM corrections in the calculation of M\ , we find 



{mff SM\ 



5rrih 



(3.12) 



In Eq.( |3.12] ), SAi^ accounts for the change in the 6-quark mass due to 0{a 2 s ) non-BLM 
corrections in JV{^ . Explicitly, 



SMf = (*» 



(2) 



r(0)r(2) 
L \ ^0 



r(l)r(l) 
^0 ^1 



L 



(0) 



L 
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L 
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L 



(i) 



(0) 



(3.13) 



To calculate d-Mf /dm&, we employ the leading order expression for A<ff , neglecting both 
perturbative and non-perturbative corrections; the impact of this derivative on the shift in 
the 6-quark mass is small. We find the following change in the value of the 6-quark pole 



mass 



drrib 



-6.6(1.6) MeV, E cut 
-6.4(4.0) MeV, E cut 



0.6 GeV; 
1.5 GeV, 



(3.14) 



where we used a s (mb) = 0.22. In brackets, the uncertainties in the mass shift related 
to numerical integration errors are indicated. Note a very strong amplification of the 
numerical integration errors when we pass from L to M moments - integration errors are 
just a few percent in the former and up to 60 percent in the latter. This implies a very 
strongly cancellation between radiative corrections in the ratio of L\ and Lq. 

Note that Eq. fEuj ) gives corrections in the pole mass scheme and that additional 
non-BLM corrections appear if the pole mass is transformed to the kinetic mass @; those 
corrections were computed in [45|. For the kinetic mass at fj, = 1 GeV, the additional shift 
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Table 4: Moment L\ in dependence of r and £. Entries at £ = are initial conditions for the 
fit; see text for details. 



is about 15 MeV, so that the total shift 

5mf n (l GeV) m 10 MeV (3.15) 

can be expected 2 . There are two ways to look at the significance of this result. We 
can compare it to the uncertainty in the 6-quark mass of about 40 — 50 MeV, typically 
obtained in fits to moments of semileptonic 6-decays [|ll|, |l2|, |l3|, ||. This comparison 



indicates that the shift shown in Eq.(3.15) is rather small. On the other hand, the error on 



the 6-quark mass in the fits is related to the fact that global fits are not very sensitive to 



2 We point out that explicit formulas that relate perturbative QCD corrections to the inclusive semilep- 
tonic b — > X c li>i decay width in the pole and kinetic schemes are given in Ref. Wx. 
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Table 5: Lepton energy moments for rrib = 4.55 GeV and m c — 1.08 GeV. 



nit, an d ni c individually; rather, the linear combination m?, — 0.6m c is restricted to about 
6 MeV. Because we estimated the shift in the 6-quark mass for the fixed value of the 
charm quark mass, a more relevant uncertainty in the 6-quark mass to compare should be 
just these 6 MeV, which is similar to our estimate of 5rrib due to non-BLM corrections 3 . 
We emphasize that the change in the 6-quark mass shown in Eq.( |3.15| ) is only an estimate 
and more careful calculation, that includes larger number of moments and simultaneous 
extraction of all heavy quark parameters, is required. 

Finally, we stress that, regardless of what uncertainty drrib should be compared to, 
Eq.( p,14j ) is remarkable since it shows that the correction to a low-energy observable due 
to two-loop non-BLM QCD effects is much smaller than the naive estimate suggests 

5m b 3 /a. 



10" d «.C F C A — ~2-10~ 2 . (3.16) 
nib V 7r / 

This feature is a consequence of a very strong cancellation between corrections to and 
L(°) ; when the ratio of the two is taken to compute Ai^ . To illustrate this point, note that 
if we set non-BLM corrections in ZA°) to zero, the shift 5rrib increases from about —7 MeV 
, as shown in Eq.( 3.14| ), to about —100 MeV. It appears therefore that high degree of 
cancellations of radiative corrections between different L-moments is crucial for claiming 
very small errors in nib, m c etc. In this respect, it is important to understand the origin 
of these cancellations since there are yet higher-order perturbative effects about which 
nothing is known at present and that, naively, are of the same order of magnitude as the 
non-BLM corrections computed in this paper. For example, although leading order BLM 
corrections C(a"/?Q _1 ) are known and resummed jnj, subleading BLM effects 0(o;!?/3q _2 ) 
are not known beyond n = 2. But, because (3q ~ 10 is large, one should expect that three- 
loop subleading BLM corrections to Lq and L\ are of the same order of magnitude as the 
two-loop non-BLM effects considered in this paper a^fio ~ a 2 s . The only way to avoid large 
shifts in the 6-quark mass is to have nearly complete cancellation between these three- and 
higher-loop corrections to Lq and L\ . The degree of such cancellation is an assumption in 
existing fits to semileptonic moments in £>-decays, as long as the origin of this cancellation 
is not clearly understood. To this end, it is interesting to give a few arguments in favor of 
non-accidental nature of these cancellations. 

For example, it is easy to see that in the limit of a very high cut on the lepton energy, 
all perturbative corrections to normalized moments vanish. Indeed, we consider the n-th 



3 We thank N. Uraltsev for emphasizing this point to us. 
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normalized moment of the lepton energy, computed in perturbation theory 
We now make a simple observation that 

^^(^ut) = (Err = {m \~ m fn n i 



independent of the strong coupling constant a s . Eq.( 3.18| ) implies perfect cancellation of 



all radiative corrections to normalized moments in that limit. Corrections to this result 
scale as 0(naJ(-E'f iax ~~ ^cut) / Ef 1 ^ ioP : ((£^ max - E cut )/rrib); they are clearly much smaller 
than the naive 0(a^) estimate of a /c-loop QCD corrections. 

Moreover, one can relax the requirement of a high lepton energy cut, by making the 
following observations. The lepton energy distribution has a peak, at Ei ~ 0.8-E™ ax . In the 
limit when this peak is infinitely narrow, normalized moments of, say, lepton energy are 
obviously protected from radiative corrections. Hence, deviations from the "no radiative- 
corrections" limit must be correlated with the broadness of the peak. To this end, note 
that the peak appears to be fairly narrow - for example, the position of the peak is only 
fifteen to twenty percent higher than the average lepton energy. We believe that results 
of explicit computations supplemented with these considerations, give a strong argument 
in favor of non-accidental nature of the observed cancellations in normalized moments for 
any value of the cut on the lepton energy and suggest that similar cancellations persist to 
higher orders in perturbative QCD. 

3.2 Decay B -> X c tv t 

We come to the discussion of the NNLO QCD corrections to semileptonic B decays into 
final states with the charm quark and the t lepton. The corresponding branching ratios 
were measured at LEP by ALEPH and OPAL collaborations [0, ||. The results read 



(2.43 ±0.32) x 1(T 2 , ALEPH, 
(2.78 ±0.54) x 10~ 2 , OPAL, 



where we added statistical and systematic errors in quadratures. We employ the ALEPH 
measurement in the following numerical computation. Using the world average for semilep- 
tonic branching ratio into the massless lepton, 

Br(B -> X c Wi) = (10.25 ± 0.25) x 10" 2 , (3.20) 

we find the ratio of the two branching fractions 

T(B^X c tv t ) , . , 

R= r ( ^«r t237(3 "' < 3 ' 21) 

The ratio 7Z can be very accurately predicted in perturbative QCD. Indeed, setting 
mi) = 4.6 GeV, m c = 1.15 GeV and m T = 1.8 GeV, we obtain the following results for 
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semileptonic decay rates 



T(b — ► ceis, 



e) 



To 

— > ctza- 



^0 (Pc Pe) 
: ^0 (PcPr) 



1 + — (-1.777) + ( — V (-1.92/?o + 3.38) 
1 + — (-1.462) + (^f (-1.82/3o + 3.16) 

7T V 7T / 



(3.22) 
(3.23) 



r 

where p c = m 2 /m 2 ,p e = 0, p T = m 2 /m 2 and To = G 2 F \ V c 6| 2 m|/(1927r 3 ). The function 
z o(Pq,Pi) re ads [||, |(| |l| 

^o(p g , P/) = (1 - 7p q - 7p 2 q + p 3 q - 7pi - 7pf + pf + p q pi (12 - 7p q - 7 pi)) 

+12^(1 - Pl f log i±^i + 12^(1 - ^ log 1+21, (3.24) 

1 — V n \--Vl 



where A(p g , pi) = 1 + p 2 + p 2 - 2p q - 2pi - 2p q pi and v q = * " , v t 



Taking the ratio, we find 



1 + p q - pi 1 + pi - p q 



npcrt = zo(pc,p T ) A + Q /oa + + /a. 

^o(PcPe) V V ^ 7 V ^ 

= Zo(Pc ' Pt) (1 + 0.0221 o(as) + 0.0044 B lm + 0.0017 non _BLM) , (3.25) 

MPcPe) ^ 

where at the last step a s = 0.22 was used. We observe that QCD effects in the 1Z ratio are 
very small. We stress that while the ratio of leading order decay rates is a rapidly changing 
function of irib,m c and m T , radiative corrections to Br(6 — * X c tv t ) and Br(6 — * X c lv{) are 
correlated, so that they cancel out in the ratio largely independent of the quark masses. We 
point out that non-perturbative corrections to the 1Z ratio were computed in |30|, [IT] and 
were found to be of the order of minus four percent. Interestingly, not only perturbative 
and non-perturbative corrections are small individually, but they also tend to cancel each 
other. Given that perturbative and non-perturbative effects in the 1Z ratio are very small, 
we fix the 6-quark mass to its value from semileptonic fits mj = 4.55 GeV and require 

Zo [ pc > pT ] = 0.237(31), (3.26) 

to determine the charm quark mass. The dependence of decay rates on quark masses at 
leading order is well-known; it can be extracted from Refs.|2^, 30, pi] ]. We obtain 

m c = 1.040(200) GeV. (3.27) 

This m c value is perfectly compatible with, but a factor of three less precise than, the 
recent result from global fits m c = 1.080(70) GeV 0. Nevertheless, the 1Z ratio seems 
to be an interesting observable since it is primarily sensitive to phase-space ratios and is 
almost independent of both perturbative and non-perturbative effects. The reduction of 
the experimental error in the 1Z ratio by a factor of three will lead to the determination of 
the charm quark mass with the precision comparable to the precision currently achieved in 
global fits. As the final remark, we point out that for central values of bottom and charm 
quark masses determined from semileptonic fits ||, mj = 4.55 GeV and m c = 1.077 GeV, 
the 1Z ratio is 0.232, in perfect agreement with the ALEPH result Eq.( 3.2l| ). 
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4. Conclusions 

In this paper, we studied the NNLO QCD corrections to semileptonic b — > X c li>i decays. We 
described the computational method that allows us to consider decays into both massless 
and massive leptons and impose arbitrary cuts on the final state particles. 

We showed that non-BLM NNLO QCD corrections to b — > X c li>i decays, with I = e, fi, 
are not very sensitive to cuts on the lepton energy, as long as the cut is below 1 GeV. 
For higher values of the lepton energy cut, the non-BLM corrections do develop E cut 
dependence, although it is not very strong. We also found that there are very efficient 
cancellations of QCD radiative corrections to normalized moments, that are used in global 
fits to semileptonic -B-decays, and that such cancellations are crucial for making the claimed 
accuracy of the fits credible. We argued that there are good reasons to believe that such 
cancellations are not accidental, and that they persist in higher orders of perturbation 
theory as well. 

We also computed QCD radiative corrections to the ratio of branching fractions of 
b — > X c tv t and b — ► X c ev e decays. It turns out that radiative corrections to this ratio 
are very small and convergence of the perturbative expansion is excellent. Since non- 
perturbative effects are also moderate, this ratio is, potentially, a good source of information 
about bottom and charm quark masses. We showed that if the charm quark mass is 
extracted directly from this ratio, the result is in good agreement with the value of the 
charm quark mass obtained from fits to semileptonic -B-decays. 
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